Power-law tails in triple system decay statistics 
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We have investigated the decay statistics of triple systems with different masses in New- 
tonian dynamics. We demonstrate that in a broad interval of mass ratios this statistics 
has good approximation by power-law tails. The power indices do not show any significant 
dependence on mass ratios. 
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[s I, though strict mathematical results about this only for a limited time interval, and for large 
kind of systems are still rather poor. enough time t power-law tails become clearly de- 



Three-body problem is an example of com- tectable |17| ]. Both these works deal with the 
bined dynamics - it is known that stable periodic equal-mass problem. This construction may be 
orbits surrounded by islands of regular dynamics too symmetric for making general conclusions 
exist in a chaotic "sea" representing the bulk of (keeping in mind loss of stability for "figure- 
phase space. For equal-mass problem (when all eight" solution" with a small mass inequality), 
three masses have the same value) there exists so it is reasonable to provide similar analysis in 
stable configuration called a "figure-eight" orbit a general situation with non-equal masses. In 
0], as well as several other orbits [lO]. It is also the present paper we have chosen several differ- 
known that the "figure-eight" orbit becomes un- ent mass ratios for investigation, starting from 
stable when relative differenses in masses reach near-equal mass system and concluding for the 
approximately 10~^ [ll], so this very interesting system with masses (0.1, 1, 10) mass units, which 
orbit is hardly important for any real astrophys- means that Schubart orbit is already unstable. 



ical problem. The Schubart orbit [13] remaines and Lagrange triangle points are still unstable 

stable in a much broader interval of masses [l^ for this mass ratio. 
On the other hand, famous triangle Lagrange 

points become stable in the situation of suffi- II. SYSTEM AND INITIAL DATA 

tiently different masses, and we can imagine that 

, r , , 1 1 1 ,1 ■■•IT In the present paper we cover wide range of 
set of stable orbits have rather nontrivial depen- 

1 mass ratios from equal-mass problem to near hi- 
dence on mass ditterences. 



The phenomenon of trajectory "stickness" to 
boundary of a regular region for the three-body 
jroblem has already been remarked by [l^ and 



10[. Keeping this in mind we can expect exis- 



erarchical system. We use Aarseth code [M] and 
start from zero- velocity initial conditions. Time 
variable is chosen so that G = 1. 

The layout of the initial conditions we used 
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tial ansatz for N{t) have been used. However, 
Orlov et al. indicated that this modelling works 



n is a space scale and 6 is some small sepa- 
ration introduced to avoid the situation when 
two bodies have totally coinciding initial coor- 
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FIG. 1: The layout of the initial conditions (coordi- 
nates). The mi body is always in the (0, 0) point, 
x-coordinate of TO2 is varying within certain range. 
As for ma, its both coordinates are varying — see 
text for details. 

dinates. Finally, the third body ms initially 
lies within a square (6, c) with b G (0, 2n) and 
c G (0 + 5, 2n + (5) . This choice cover rather wide 
area in the initial condition space for our results 
to be representent, on the other hand, different 
permutations of non-equal masses results in dif- 
ferent initial condition sets allowing us to study 
possible dependence of our results upon initial 
conditions chosen. 

Using the described above layout we evenly 
distribute about 4 millions initial configurations 
and calculate their evolution until the decay (for 
exact meaning of the decay in Aarseth code see, 



for example, 




FIG. 2: The typical plot of N{t) in the logarithmic 
scale. 

III. DECAY SPECTRA 

A typical plot of the differential distribution 
N{t) in the logarithmic scale is shown in Fig. 
2. Power-law tail following a steeper function 
after t reaches the value of about one thousand 
is clearly seen. The distribution of escaping time 
for small t is clearly different from power-law, 
and to obtain correct values of power index in 
the modelling 

N{t) = const X r"/ 

we should exclude this range of t from the anal- 
ysis. As an illustration we plot in Fig. 3 power 
indices calculated in the time range from tm 
to Tf where the maximal time of integration 
Tf = 7000. We can see that the power index 
typically does not change if tm is bigger than 
2000. In what follows we will ignore first 2000 
time units in calculation of power index 7. In 
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the table below we present minimal and maximal 
values of 7 for several mass combinations calcu- 
lated for different sets of initial conditions. Error 
bars, also presented in the table indicate that in 
some situations the difference in power indices 
calculated for the same system using different 
sets of initial conditions can be real (though not 
so big). We can see also that despite mass ratio 
of the systems under investigation cover a broad 
interval from near-equal to near-hierarchical sys- 
tems, the value of 7 does not change significantly 
with the slightly developed tendency of growing 
7 with the increasing hierarchicity of the triple 
system. We should, however, note that our val- 
ues of 7 for equal mass system is bigger than 
the value found in the paper [17]. Can different 
sets of initial conditions be the only cause of this 
difference remaines unclear and requires further 
investigations. 



row interval (from 1.76 to 2.19). 

TABLE I; Maximal and minimum values for 
the slope 7, calculated at tin ~ 3000 for dif- 
ferent mass combinations. First column gives 
absolute masses, second - masses, normalized 
to mi + m2 + fTi3 = 1, third - the slope 7 with 
an error, minimal (upper row) and maximal 
(lower row) values. 





Masses 


7 ± ^7 


(0.1, 1, 10) 


(0.009, 0.09, 0.901) 


1.762 0.011 
1.863 0.016 


(1, 3, 10) 


(0.071, 0.214, 0.715) 


1.845 0.013 
1.862 0.012 


(1, 3, 5) 


(0.111, 0.333, 0.556) 


-1.828 0.013 
1.848 0.014 


(1, 2, 3) 


(0.167, 0.333, 0.5) 


1.857 0.014 
1.905 0.015 


(0.9, 1, 1.1) 


(0.3, 1/3, 11/30) 


2.149 0.009 
2.181 0.010 



(1, 1, 1) (1/3, 1/3, 1/3) 2.125 0.014 

2.178 0.012 



IV. CONCLUSIONS 

We studied the escaping rate statistics in 
a general three-body problem with non-equal 
masses. We argue that for large enough time 
t this statistics can be modelled by power-law 
functions with a good accuracy. For all cases 
studied (with normalized mass ratios from to) 
the power index is located within a rather nar- 
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FIG. 3: An examples of dependence of the slope 7 on the cutoff time Black lines corresponds to (1, 2, 
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